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Abstract. Over the last years, lattice calculations in pure Yang-Mills gauge theory seem to have come more or less to a 
consensus. The ghost propagator is not enhanced and the gluon propagator is positivity violating, infrared suppressed and 
non-vanishing at zero momentum. From an analytical point of view, several groups are agreeing with these results. Among 
them, the refined Gribov-Zwanziger (RGZ) framework also accommodates for these results. The question which rises next 
is, if our models hold the right form for the propagators, how to extract information on the real physical observables, i.e. the 
glueballs? How do the operators which represent glueballs look like? We review the current status of this matter within the 
RGZ framework. 
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INTRODUCTION ON THE GZ 
FORMALISM 

As is well known, the Yang-Mills action needs to be 
gauge fixed in order to define the path integral. One way 
to do this is trough the gauge fixing proposed by Faddeev 
and Popov. However, as is shown by Gribov, this gauge 
fixing suffers from Gribov copies [1]. If we take e.g. the 
Landau gauge, d^A^ = 0, it is very easy to show that 
there exist gauge equivalent fields A'^ also fulfilling the 
Landau gauge d^A'^ when having zero modes of the 
Faddeev-Popov operator 

J(* = -d ll D ll = ~d ll {d ll d" h ^gf' bc A^). 

Therefore, Gribov proposed to restrict the region of inte- 
gration to Q. 

Q = {A tt |^A /i = 0,^>0)} , 

the region which is free of zero modes of the Faddeev- 
Popov operator ^ ab . We should mention however, that 
there are still Gribov copies inside the Gribov region 
[2]. In 1989, Zwanziger implemented this restriction to 
the Gribov region to all orders: the Gribov-Zwanziger 
action was born [3]. Let us here immediately present the 
localized version of this action, namely 

Sqz = Sfp + So + Sy , 



with Spp the usual Faddeev-Popov action, 

S FP = \j d 4 F« v F« v + J d 4 x (b a d^+^Dfc b ) , 

(1) 

and Sq and Sy given by 

S = J d 4 x {l?™J( ah ($ - (Sgf.*"*®*) , 

Sy = -fg J d 4 x (r 6c (<+^ c )A£ + - g ^ 2 - ! ) A • 

The fields (<p^ r ,<p^ c ) are a pair of complex conju- 
gate bosonic fields, while (7o^ c , G)!f) are anticommuting 
fields. 7 is not free, but fixed by a horizon condition, 
{gfdJfytff + ffi) = 2yH{N 2 - 1). Looking at this 
horizon condition, we see that it is equivalent with a di- 
mension two condensate. Therefore, we could investigate 
other dimension two condensates, namely (<p<p — Wco) 
and (A 2 ). Including these condensates gives rise to the 
refined GZ (RGZ) action [4]. We have shown that this 
can be done in a renormalizable way. We stress that the 
masses related to the condensates respectively M 2 and 
m 2 are dynamically generated. 

PROPAGATORS 

Let us now look at the ghost and the gluon propaga- 
tor. Firstly, it is sure that perturbation theory fails in the 



infrared region as perturbatively, the gluon propagator 
blows up, while all lattice calculations show that there 
is infrared suppression. Secondly, the Gribov-Zwanziger 
predicts a gluon propagator which is infrared suppressed, 
and therefore clearly goes beyond perturbation theory. 
Unfortunately, at zero momentum, this gluon propagator 
vanishes, which contradicts the current numerical simu- 
lations that this propagator is non-vanishing at zero mo- 
mentum [5, 6]. In addition, the ghost propagator shows 
infrared enhancement, which is also no longer observed 
in the infrared. Therefore, something is still missing in 
the GZ framework. Thirdly, let us investigate the refined 
GZ action. In this case, the tree level gluon propagator 
behaves like 



representation: 



D(p 2 ) = 



p 2 +M 2 



(M 2 - 



»2\ „2 



p 2 + 2g 2 Nf+M 2 m 2 



(2) 



This propagator is clearly infrared suppressed, also vio- 
lating positivity and is no longer zero at zero momentum. 
Also the ghost propagator is no longer enhanced as an 
effect of the condensates. This agrees again qualitatively 
with the latest lattice results. 

In [7], we have investigated our specific form of the 
gluon propagator on the lattice. Two main conclusions 
could be drawn: 1) If m 2 = 0, the fits are of low qual- 
ity, therefore, the condensate (A 2 ) ~ m 2 is indispensable; 
2) Nice fits for the form (2) were found when includ- 
ing the condensate (A 2 ). Let us recall some estimates: 



2.14±0.13GeV 2 , m z 



1.78 ±0.14 GeV 2 and 



D(0) = 8.2±0.5 GeV 2 . From the estimate of m 2 , we find 
that (<? 2 ^ 2 ) M= ioGeV ~ 3GeV 2 which is in the same ball- 
park as other, rather independent, approaches [8, 9]. 
In conclusion, we could state that the RGZ framework 
provides a possible explanation of the behavior of the 
gluon and the ghost propagator. In fact, one can show that 
more condensates can be present, and alternative RGZ 
frameworks exist [10, 11]. 

Let us also stress that in 2d, the ghost propagator is still 
enhanced, while the gluon propagator is zero at zero mo- 
mentum, [12]. We have also provided evidence that in 
2d, the GZ formalism still holds the right results as re- 
finement is impossible in 2d [13]. 



THE QUEST FOR PHYSICAL 
OPERATORS 

Now that the propagators match with the lattice, we can 
wonder whether the RGZ also holds information about 
the particles of (quenched) QCD? The idea to proceed 
is the following: we want to find an operator G, so that 
the correlator (0{k)0{—k)) can be put into a spectral 



/■CO 

(0(k)0(-k)) = / dxp(x) 



x+k 2 



JTdr p(r) ^ so that when 



&5L.I.. we find 



We then introduce F(z) 

going to Minkowski space, i.e. k\ ud -r -, Mink 
a discontinuity along the positive real axis. In addition 
we want p(x) to be positive for x > Xq in order to give 
a particle interpretation to the correlator 1 and obviously 
0{k) has to be renormalizable. 

We have worked out two different angles to attack this 
problem. In the first approach, we have introduced i- 
particles [14] 



Jf A ^2^ r "(< + ^)' 



to make the quadratic part of the action diagonal. With 
the /-field strengths, X£ v = d^X" — d v X£ and r]^ v = 
dfiT]y — dyX)^, we have proposed the following operator 

The corresponding correlator (0^ l \k)0^(— k)) can be 
put in a spectral representation: 

r°° 1 
{ 0^\k)O^H-k)) = j n2 dxp(x)—,, 

whereby p(x) > 0. Unfortunately, the operator is 
not renormalizable. 

Let us therefore also discuss the second angle. In [15], 
we have investigated the renormalization of the typical 
scalar glueball operator F 2 with the GZ action. Let us 
start by noting that the Faddeev-Popov action Sfp is 
invariant under the BRST transformation s: sSw = 0, 

sA; = -( Dfl c) a , sc a =\gf bc c h c c , 



sc a = b a , 



sb a = . 



This BRST invariance lies at the origin of the Slavnov- 
Taylor identity, which allows us to prove the renormal- 
izability. Moreover, the BRST charge allows us to define 
the sub-space of the physical states and to establish the 
unitarity of the S matrix. Unfortunately, the GZ action 
is no longer invariant under the BRST transformation s. 
The new fields transform as 



03; 



«9f = 0, «9? = pf, sff = 0, (3) 



p is proportional to the cross section, and thus has to be positive. 



and thus the breaking is proportional to y, sSgz = 
s (Sym + S g f + S + S 7 ) = s (S r ) ~ y 2 ^ 0. Despite this 
breaking, the GZ action is still renormalizable, due to 
a rich set of Ward identities. Moreover, only two renor- 
malization constants are needed, which is the same as in 
the Yang-Mills theory. Let us also mention that by in- 
troducing extra fields, a symmetry can be restored again 
[16]. With the breaking in mind, we can renormalize F 2 
within the GZ framework as done in [15]. This renormal- 
ization is however far from trivial due to the breaking of 
the BRST symmetry. We find that F 2 mixes with the fol- 
lowing operator: 

S = s(. . .) + fDf +d(N 2 - 1) 7 . (4) 

Moreover, we have constructed a renormalization group 
invariant given by 



(2) = ^l F 2_ 2Yc£ 

8 



(5) 



Therefore, we would propose the following correlator 



{2) {x)O i2) {y) 



■2y c < 



to 



^F 2 -2y L , 



(y) 



* (F 2 ( X )F 2 (y)) . 



(6) 



Notice that the breaking of the BRST symmetry is the 
reason that this correlator contains extra terms beside 
(F 2 (x)F 2 (y)Y In ordinary Yang-Mills theory, a similar 
correlator can be found, but it reduces the (F 2 (x)F 2 (y)} 
due to the presence of the BRST symmetry. Unfortu- 
nately, with this correlator, we are unable to find a good 
spectral representation. This was already found in [3]: 

J d 4 xe- ikx (F 2 (x)F 2 {0)) = G phys (k 2 ) + G unphys (k 2 ) , 

whereby the unphysical part, G unphys (£ 2 ), displays cuts 
along the imaginary axes beginning at the unphysical 
values k 2 = ±4iy 2 and the physical part, G phys (k 2 ), has 
a cut beginning at the physical threshold k 2 = —2y 2 . 

For further clarification, let us show the relation be- 
tween M 1 ^ and M 2 \ At lowest order, we have that 



OP) 



1 a n a 
"-fiv'lfiv 



and thus it are in fact the two last terms which cause 
the bad spectral representation in (F 2 (x)F 2 (y)), and the 
question remains how we can find a renormalizable op- 
erator with a good spectral representation. 



We can already mention that so far, we have only investi- 
gated spectral representations within the GZ framework. 
The hope is that within the RGZ framework, we can still 
find an operator with the required properties as some as- 
pects change of the analysis. Also, the study of glueballs 
looks very promising within the (R)GZ framework, see, 
for example, the recent results in [17]. This question of 
finding descent operators is not only relevant within the 
Gribov-Zwanziger context, but for all people involved in 
infrared propagators QCD. How to find good spectral op- 
erators starting from unphysical operators? 
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